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ABSTRACT 

The interstellar medium is not necessary fully ionized. When a shock propagates into the partially 
ionized medium, some neutral particles leak into the upstream region from the downstream region. In 
this paper, we investigate how leakage neutral particles affect the upstream structure of the shock and 
particle accelerations. By using four fluid approximations (upstream ions, upstream neutral particles, 
leakage neutral particles and pickup ions) , we provide analytical solutions of the precursor structure 
due to leakage neutral particles. It is shown that the upstream flow is decelerated in the precursor 
and the shock compression ratio becomes smaller than that of the strong shock limit, but the total 
compression ratio is the same as that of the strong shock limit. Even if leakage of neutral particles is 
small (a few percents of upstream particles), this smaller compression ratio of the shock can explain 
steep gamma-ray spectra from young supernova remnants. 

Subject headings: acceleration of particles — cosmic rays — plasmas — shock waves — supernova 
remnants 



1. INTRODUCTION 

f Supernova remnants (SNRs) are thought to be the 
origin of Galactic cosmic rays (CRs). The most popu- 
lar SNR acceleration mechanism is diffusive shock ac- 
celeration (PSA) LA xford et al. 1977; Krv mskvl 119771 : 
[Belil[l97a iBlandford fc Ostrikeri Il978l) . In fact, Fermi 
and AGILE obs erved GeV gamin a rays due to CRs 
from SNRs fe.g. lAbdo et all T200l [2015 iTavani et al.l 
[20T0I: iGTuliani et all 120111: IGiordano et all 120121) . How- 
ever, gamma-ray spectra from SNRs are steeper than 
that expected from the standard DSA theory. The 
steep spectra can be interpreted as effects of a n en- 
ergy dependent escape jP tuskin & Zirakas hvilil 120051 : 
lOhira et al.ll2010l:lCaprioli~t al...2Q10. ) and diffusion (e.g. 
lOhira et alllMll) for middle-aged SNRs (- 10* yrs old). 
For young SNRs (< 10^ yrs old), some ideas to ex- 
plain the steep spectra have b een p r oposed (Kirk et al.l 
1996; Zirakashvili & Ptusk inI 120091: lOhira et al.l 120091: 
Ohira fc Takaharai201Qi: ,Bell"et al.) 120121) but it is still 



an op en issue. " lOhira et al.l (|2009f ) nOhira fc Takaharal 
(|2010D proposed importance of neutral particles to parti- 
cle accelerations and shock structures. A neutral fraction 
of the interstellar m edium around SNRs is often found 
to be order of unity ([Ghavamian et al.ll200a I2002D . 

Interactions between neutral particles and plasma have 
observed directly in the solar wind. The re are two pecu- 
liar p articles, energetic neutral atoms (iMcComas et al.l 
120091) and pickup ions (Gl oeckler et al.l ll993') in the so- 
lar wind. Although origins have not been completely 
understood, an attractive idea has been proposed. Neu- 
tral particles penetrate into the inner solar system from 
the surrounding local interstellar medium. The neutral 
particles have large velocity in the rest frame of the so- 
lar wind. After they are ionized by charge exchange and 
photo ionization, they are picked up by the magnetic field 
of the solar wind, so that they become pickup ions with 
large velocity dispersion in the solar wind. After passing 
over the termination shock, some pickup ions become 
neutral atoms by charge exchange. The neutral atoms 



can propagate toward the sun and are observed as ener- 
getic neutral atoms. 

Applying these pictures to SNR shocks propagating 
into a partially ionized medium, we expect leakage of 
neutral particles from the downstream region to the up- 
stream region. Furthermore, the leakage neutral particles 
become pickup ions by coUisional ionization or charge 
exchange in the upstream region. If leakage is signif- 
icant, the upstream flow is decelerated and heated by 
the pickup ions. As the result, a precursor is produced 
by lea kage neutral particles. Very recently, iBlasi et all 
(|2012D showed that neutral particles really leak into the 
upstream region from the downstream region by solving 
the Vlasov equation of neutral particles. They showed 
the formation of the precursor due to the neutral par- 
ticles and provided the precise velocity distribution of 
neutral particles by assuming that only the ion distribu- 
tion is a Maxwellian. 

In this paper, we investigate the precursor structure 
by a different approach which is a four fluid approxi- 
mation. We consider upstream ions, upstream neutral 
particles, leakage neutral particles and pickup ions, re- 
spectively. The velocity distribution of ions in the down- 
stream region has not been understood for SNR shocks, 
especially, in the partially ionized medium. Therefore, 
we treat physical values concerning to the leakage neu- 
tral particles as free parameters. Then, we obtain some 
analytical solutions of the precursor structure. Our re- 
sults in th i s pap er are qualitatively similar to that of 
IBlasi et al.l ([2012I ). 

We first provide some lengthscales for collisional ioniza- 
tion and charge exchange in Section 2. We then provide 
four fluid models to describe the precursor structure in 
Section 3, and provide two approximate solutions in Sec- 
tion 3.1 and 3.2. Section 4 is devoted to the discussion. 

2. RELEVANT LENGTHSCALES 

In this section, we briefly summarize relevant length- 
scales in the precursor due to leakage neutral particles. 
We here consider charge exchange and collisional ioniza- 
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tion as interactions between ions and neutral particles. 
Although collisional ionizations with electrons could be 
important, it depends on an el ectron temp erature which 
has not been understood y et (jOhira iLTakahara! I2007L 
120081 iRakowski et all 120081 ). We here do not take into 
account ionization due to electrons. Cross sections of 
charge exchange and colli sional ionization depen d on a 
relative velocity, Uroi fe.g. lHeng fc McCravll2007[ ). CoUi- 
sional ionization is dominant for u^ei ^ 3000 km s~^ and 
its cross section of hydrogen is typically di ^ 10~^^ cm^. 
On the other hand, charge exchange is dominant for 
Wrei ^ 3000 km s^^ and its cross section of hydrogen 
is about crc.E. ~ 10^^^ cm^. 

In the shock rest frame, leakage neutral particles typ- 
ically have a velocity comparable to the shock velocity, 
wieak ~ —Ush- Then, the relative velocity between up- 
stream particles and leakage neutral particles is about 
Wroi,icak ~ 2Micak- The shock velocity is typically Ush ~ 
3000 km s""'^ for young SNRs, so that collisional ioniza- 
tion with upstream ions and upstream neutral particles 
is the dominant process for leakage neutral particles. 
Therefore, the precursor scale in the shock rest frame, 
Lpre, is given by 

Lpic 5 X 10^^ cm 

rcl,lcak 

05 



( ^ 



1 cm^-^ 



(1) 



where ritot is the total number density of the upstream 
ions and neutral particles. 

Upstream neutral particles are ionized by collision with 
leakage neutral particles and pickup ions. The interac- 
tion length in the shock rest frame, Li, up, is given by 



Li,up-5 X 10^6 cm 

^sh / ^rcljcak 

(15 
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where riioak is the number density of leakage neutral par- 
ticles and we assume that the number density of leakage 
neutral particles is smaller than upstream total number 
density. Therefore, as long as Ush/'^icak > "icak/n-tot 
(ii,up > Lprc), the ionization of upstream neutral parti- 
cles by leakage neutral particles is negligible in the pre- 
cursor. 

Upstream neutral particles interact with not only leak- 
age neutral particles but also upstream ions. The rel- 
ative velocity between upstream neutral particles and 
upstream ions, Wroi.up, would becomes larger than their 
thermal velocity but smaller than the shock velocity, 
Ush ^ 3000 km s~i. Therefore, charge exchange is dom- 
inant and the interaction length between upstream ions 
and upstream neutral particles in the shock rest frame, 
Lc.E.,up, is given by 



ic.E.,up~2 X 10^^ cm 

'^sh/ '^rcl,up 

100 



\ / crc.E. \-i / rtion \ - 
J V 10-15 cm2/ Vo.5 cm-3/ 



(3) 



where riion is the number density of upstream ions. If 
leakage is small, the upstream plasma flow does not 
change significantly and ic.E.,up > ipre- Then, up- 
stream neutral particles rarely interact with upstream 



ions in the precursor, that is, the fiow velocity of up- 
stream neutral particles does not change. We consider 
this decoupling case in Section 13.11 If leakage is large, 
the upstream plasma flow is significantly decelerated and 
Urci.up ~ and -/jc.E.,up < LpiQ. Then, upstream neu- 
tral particles interact many times with upstream ions in 
the precursor, so that the flow velocity of upstream neu- 
tral particles is almost the same as that of upstream ions. 
We consider this tight-coupling case in Section [3?2] 

According to the DSA theory, accelerated particles dif- 
fuse into the upstream region. The diffusion length is 
given by 



SeBush 
■10i5 cm 



3000 km s- 



B 



100 



E 



1 TeV 



(4) 



where rjg, B and E are the gyrofactor, the magnetic field 
and the energy of accelerated particles. Therefore, parti- 
cles are accelerated in the precursor due to leakage neu- 
tral particles, that is, leakage of neutral particles is im- 
portant for the particle acceleration. The diffusion length 
of pickup ions is much smaller than the precursor scale. 



^prc- 



Hence, we can neglect diffusion of pickup ions. 



3. FOUR FLUID MODEL IN A NEUTRAL PARTICLE 
PRECURSOR 

In this section, we calculate the steady-state precursor 
structure due to leakage neutral particles. There are cold 
ions and neutral particles in the upstream region. In ad- 
dition, we consider leakage neutral particles and pickup 
ions originated from leakage neutral particles. We here 
adopt a four fluids model to describe the precursor. Con- 
tinuity equations of the steady state are given by 



dx 



(nnM„) = -Qi, 



dx 



("-IcakUL 



cak } 



dx 
d_ 

dx 



(^ion^if 

(npuiWi, 



— Qpui , 

' Qion ; 

' Qpui , 



(5) 
(6) 

(7) 
(8) 



where the subscripts 'n', 'leak', 'ion' and 'PUF denote 
upstream neutral particles, leakage neutral particles, up- 
stream ions and pickup ions, respectively and x is the co- 
ordinate of the direction along the shock normal {x = 
and X — —oo are the position of the shock and the 
far upstream region, respectively.) We assume that 
the fluid velocity of pickup ions is the same as that of 
upstream ions because of electromagnetic interactions, 
upui = Wion- Qion and Qpui are source terms and given 
by 

Qion = ("-leak + "-PUl)"-nO'iWrol Jcak , (9) 
QPUI = ("n + n-PUI + "ion) ^leak*^i^rel,leak ? 

(10) 

where we approximate all relative velocities as a constant 
(see Equation p^ ). The first and second terms of Qion 
are due to the collisional ionization of upstream neutral 
particles with leakage neutral particles and pickup ions. 
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Fig. 1. — Densities and velocity profiles for the decoupling ap- 
proximation. The green, blue, cyan and magenta lines show densi- 
ties of upstream neutral particles, leakage neutral particles, up- 
stream ions and pickup ions, respectively. The red and black 
lines show the ion flow velocity and the Mach number, respec- 
tively. Boundary conditions are Mg = 100, »in,o/{"ion,o + "-n.o) = 

0.2,nieak,sh/("ion,0 + '^n.o) = 0.05 and Mleak,sh/«0 = -0.5. 

respectively. The first, second and third terms of Qpui 
are due to the collisional ionization of leakage neutral 
particles with upstream neutral particles, pickup ions 
and upstream ions, respectively. Moreover, we here ne- 
glect charge exchange of leakage neutral particles. This 
approximation is valid for shock velocities larger than 
about 1500 km s~^. For the decoupling (Lc.E.,up > ipro) 
and tight-coupling (Lc.E.,up < -^prc) limits, the velocity 
difference between upstream ions and upstream neutral 
particles is small. Therefore, charge exchange between 
upstream ions and upstream neutral particles does not 
change anything. In the next subsection, we solve Equa- 
tions ([SI)-® by using the momentum and energy conser- 
vations. The boundary conditions are as follows: 



(11) 

(12) 

npui(-oo)=0 , (13) 

"-lcak(O) =nicak,sh ,'"loak(0) = Wloak.sh , (14) 

where the subscripts '0' and 'sh' represent quantities 
at the far upstream region and at the shock, respec- 
tively. Note that Mieak,sh is negative. All quantities are 
nor malized by uq and nion,o + "-n.o in Sections 13.11 and 
13.21 Hence, free parameters are three, n.„_o/(?iion.o + 
'^n,o),'«ioak,sii/(nion,o + ^^n,o) and uioak.sh/wo- Wc treat 
leakage neutral particles as a cold fluid in this paper, 
that is, we neglect the velocity dispersion of leakage neu- 
tral particles. The relative velocity, Urd, is given by 

Urel,loak = ""O " ""loak.sh • (15) 

3.1. Decoupling Approximation (ic.E.,up > -^proj 

In this subsection, we solve equations of the four fluid 
system using the decoupling approximation correspond- 
ing to small leakage of neutral particles. We here assume 
that fluid velocities of upstream neutral particle and leak- 
age neutral particles do not change from that at the far 
upstream region, uq and that at the shock, Micak,sh, re- 
spectively. Then, continuity equations are given by 



Wfeak.sh^ ("leak) — — QpUI 



(16) 
(17) 



Fig. 2. — The fluid velocity of upstream ions at the shock 
(solid lines) and the spectral index of accelerated particles (dashed 
lines) for the decoupling approximation. The red, green, blue 
and black lines show cases of (nn^o/(nion,o + «n,o), "loak,sh/"o) = 
(0.1, -0.5), (0.1, -0.4), (0.1, -0.6)' and (0.5, -0.5), respectively. 

~] ("ion'^ion) — Qion , (1^) 

ax 

(npuiMion) = Qpui • (19) 

da; 

The momentum conservation law of upstream ions and 
pickup ions is given by 



dx 



{m (riion + npui) + Pjon + Ppvi} 

= mUoQion + ™Uloak,sh(9pUI , 



(20) 

where m, Pion and Ppui are the particle mass, pressures 
of upstream ions and pickup ions, respectively. The en- 
ergy conservation law of upstream ions and pickup ions 
is given by 



1 

dx 1 2 



-TO(ni, 



"PUljUk 



7 



7-1 



'^ion 



Icak.sh 

Qpui)(2;) 



where 7 is the adiabatic index. To make the ex- 
pression simple, hereafter velocities, densities, pressures 
and spatial coordinate are normalized by uq, (T^-ion.o + 

nn,o),m{nion,o + nnfl)ul and (nion,o + ?^n,o)"^c^^^ respec- 
tively. Normalized quantities are denoted with a bar. 
From Equations ([T6| - (|2ip . one can make four conserved 
quantities. Hence, Equations (ITB|) - (I^ reduce to follow- 
ing two equations: 



da: 



-UrGl,loak"nn-lcak 1 



d _ Urol, leak - /-I 

-T^»^lcak = »^lcak "-n H 

dx Uioak,sh V 

where be expressed by 



^IcakjSh 
^ion 

71y^ ^Icak^lcak.sh 
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ion 



_B + - 4AC 

ion- ^ , 

^ = (7 + 1) (1 - fin - nioak^^loak.sh) 



(22) 

23) 

(24) 



B = 27 1 



C = (7-l) 1 + 



^^j2 ^loakWjcak.sh 
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1)M2 
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4 



Yutaka Ohira 



where Mq is the Mach number at the far upstream region. 
From Equations (l22]) - (p4)) . one can obtain solutions of 
n^{x) , n\ciik{x) and Uion(a;) by numerical computations. 
By assuming that the brackets terms of Equations (|22p 
and (j23p are constant, one can obtain analytical approxi- 
mations which are valid for small fiicak.sh- Using the solu- 
tions of nioak and Mion, the other quantities, nion, "-pui 
and Pion + Ppvi can be expressed by 



_ _ I - Tin 

^ion ^ _ : 
^ion 

_ ^lcak,sh _ 

^PUI — ^Icak , 

^ion 

- 1 _ _ _2 

-PPUI — IH 772 - "n - "IcakWicak.sh 



(25) 
(26) 



- (1 - fin - nioakUloak,sh) ""ion • (27) 

Moreover, the evolution of the Mach number, M{x), can 
be expressed by 



M = uu 



72 




(•28) 



Figure [T] shows numerical solutions to Equations 
- (1251), where Mq = 100,7 = 5/3, fin.o = 0.2,nicak,sh = 
0.05, and uicak.sh — —0.5. For the same input param- 
eters, the analytical approximations of Equations (j22l) 
and ([23l) give an about 0.1 percent accuracy. There is 
no solution with M < 1. The flow velocity of ions and 
pickup ions, Uion, is slightly decelerated by small leakage 
of neutral particles from the downstream region. How- 
ever, the Mach number becomes small significantly be- 
cause the pressure of pickup ions is lar ge. The se features 
are qualitatively the same as results of (jBlasi et al. 2012). 
The number density of upstream neutral particles, is 
almost constant. As already mentioned in Section 2, ion- 
ization of upstream neutral is negligible in the precursor 
for the small nioak,sh- Therefore, the number density of 
upstream neutral particles, fin, can be approximated by 
constant. Then, one can obtain analytical solutions at 
the shock. 

Because we have all quantities at the shock, we can cal- 
culate the shock jump condition. The collisionless shock 
is formed by only the plasma because the dissipation 
length of the plasma is much small than the interaction 
length of neutral par ticles (Chevalier & Raymond 1978; 
iChevalier et al.l [19801 ). It should be noted that we have 
to take into account sink terms due to leakage of neutral 
particles when we derive the Rankin-Hugoniot relations. 
Even though leakage neutral particles are not ions, the 
origin is hot ions in the downstream region. By using 
the Rankin-Hugoniot relations between the far upstream 
region and the down stream region, one can easily ob- 
tain the shock jump condition. This is because the net 
source of ions is only upstream cold neutral particles in 
between the far upstream region and the downstream re- 
gion. Therefore, the total compression ratio between the 
far upstream region and the downstream region is given 




0.01 



-0.5 

X / ("ion.O+nn.o)"! 



Fig. 3. — The same as Figure [T] but for the tight-decouphng 
approximation. Boundary conditions are Mq = 100, n„_Q/(nion.o + 
nn,o) = 0.9,nieak,sh/("ion,o + "n,o) = 0.3 and ■Ulcak,sh/'fO = -0.5. 



by 



r-tot 



(7 + 1) 



(7 - 1) + 2Mo 



(29) 



where rtot = 4 for 7 = 5/3. The partial compression 
ratio between the far upstream region and just behind 
the shock — —e) is Uj^^, so that the compression ratio 
of the shock is given by 



rsh 



(7 + I) 



(7-l) + 2Mo 



i(i = 0) 



(30) 



One can obtain the analytical approximations of the 
shock compression ratio, rsh, because of n„ w nn,o- Fig- 
ure [2] shows the analytical approximations of the fluid 
velocity of upstream ions at the shock, Uion{x = 0), 
and the spectral index of accelerated particles, s = 
(rsh + 2)/(rsh — 1)- Even though the number density 
of leakage neutral particles is small, nioak.sh ^5 the 
the spectral i ndex beconi e s larg er than 2. As already 
mentioned by iBlasi et "aD (|2012j ). this can explain the 
observed gamma-ray spectra slightly steeper than the 
simplest prediction of DSA. Effects of leakage neutral 
particles becomes large when the ionization fraction is 
small and the flux of leakage neutral particles is large. 
Note that the decoupling approximation becomes some- 
what bad for fiioak,sh ^0.1 because the relative velocity 
between upstream neutral particles and upstream ions 
becomes large (see Section 2). 

3.2. Tight- Coupling Approximation (lc.E.,up < -^prc/ 

In this subsection, we solve equations of the four fluid 
system using the tight-coupling approximation corre- 
sponding to large leakage of neutral particles. We here 
assume that the fluid velocity of leakage neutral parti- 
cles does not change from that at the shock, wieak.sh- In 
addition, we assume that the fluid velocity of upstream 
neutral particles is the same as that of upstream ions be- 
cause of the tight coupling, Un = Uion- Then, continuity 
equations are given by 



ax 



"lGak,sh^ (,«lcakj 



-Qpui 



lon^ion 



(31) 
(32) 
(33) 
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Fig. 4. — The same as Figure [2] but for the tight-decoupUng 
approximation. The red, green and blue Unes show cases of 
jileak,sh/fO = "0.5,-0.4 and —0.6, respectively. 



da; 



(npuiWion) = Q: 



PUI 



(34) 



The momentum conservation law of upstream neutral 
particles, upstream ions and pickup ions is given by 



d 

da; 



{m (n,! + nion + n-pui) wfon + + P\on + Ppvi] 

= 77micak,shQpUI (35) 

where Pn is the pressure of upstream neutral particles. 
The energy conservation law of upstream neutral parti- 
cles, upstream ions and pickup ions is given by 



1 

dx 1 2 



m(nii + riion + npui)uf 



3 

ion 



7-1 



2"^^loak,sii'3PUI 



^ion {Pn + -Pion + ^PUl) 



(36) 



To make the expression simple, hereafter velocities, den- 
sities, pressures and spatial coordinate are normalized in 
the same way as the previous subsection. From Equa- 
tions ([3T |) - ([351) . one can obtain the following equations: 



dx 



nn — — '"rol,loak"n?^loak 1 



d _ _ Ui.el,leak _ 
~r3^1cak — ~ ^Icak 
dx Uioak,sh 



Wlcak,sh 
^ion 

1 ^lcak^lcak,sh 



Us, 



where be expressed by 



_B + - AAC 

ion- 1^ , 

A = (7 + 1) (1 - nicakWlcak,sii) , 

1 - -2 

—772 '^Icak^lcak.sli 



(37) 
(38) 

(39) 



B = 27 1 



C=(7-l) 1 + 



(7-l)Af2 



"-lcakWieak,sii 



From Equations (I571) - ([M)) . one can obtain solutions of 
nn(a;), fiicak(a;) and Uion(a;) by numerical computations. 
By assuming that the brackets terms of Equations (|37p 
and (1381) are constant, one can obtain analytical approx- 
imations. Using the solutions of nn,riieak and Wion, the 



other quantities, riion, "-pui and Pn + -Pion + -Ppui can be 
expressed by 



1 

^ion 

_ ^leak,sh _ 

npui = z ^^Icak , 

^ion 

1 _ _2 

Pn + Pion + PpUI = 1 + , ^2 ~ "-IcakMleak sii 



(40) 
(41) 



7M2 

- (1 - ?lloakUloak,sh) Wion ■ (42) 

Moreover, the evolution of the Mach number, M(a;), can 
be expressed by 



75 




m {tin + n-ion + »PUl) 
7 (Pn + Pion + PpUl) 



^ilcakWioak.sli 



^lcak'^lcak,sli) '^io 



- 1 



(43) 



Figure |3] shows numerical solutions to Equations (P7)) - 
where Mq = 100,7 = 5/3,n„,o = 0.9,nioak,sh = 0.3, 
wioak,sii = —0.5. For the same input parameters, the an- 
alytical approximations of Equations ([57)1 and ([55)) give 
an about 5 percent accuracy. There is no solution with 
M < 1. Unlike the decoupling approximation, the flow 
velocity of upstream neutral particles, upstream ions and 
pickup ions, Uion, is significantly decelerated because of 
large leakage of neutral particles from the downstream 
region. The number density of upstream neutral parti- 
cles, fin, and upstream ions, nion, become large because of 
compression. Moreover, the increase rate of the number 
density of upstream ions is larger than that of upstream 
neutral particles because some of upstream neutral par- 
ticles are ionized in the precursor. Unlike the decoupling 
approximation, solutions of the upstream ions velocity 
and the Mach number, Equations (p9| and (j43|) . do not 
explicitly contain the density of upstream neutral parti- 
cles. Therefore, one can obtain analytical solutions at the 
shock. From Equation (I5(I1) . one can obtain the analyti- 
cal solution of the shock compression ratio, r^h- Figure|4] 
shows the analytical solutions of the fluid velocity of up- 
stream ions at the shock, Mion(a; = 0), and the spectral 
index of accelerated particles, s — (rgh + 2)/(rsh — 1)- 
When the number density of leakage neutral particles is 
"■icak,sh ^ 0.1, the spectral index becq r nes la rger than 
2. As already mentioned bv iBlasi et al.l (|2012[ ). this can 
explain the observed gamma-ray spectra slightly steeper 
than the simplest prediction of DSA. Note that the tight- 
coupling approximation becomes bad for nicak,sh <C 0.1 
because the fluid velocity of upstream ions does not sig- 
nificantly change and the relative velocity between up- 
stream ions and upstream neutral particles is still small 
(see Section 2). 

4. DISCUSSION AND SUMMARY 

We here discuss some important effects neglected 
above. When leakage neutral particles are ion- 
ized, their velocity distribution in the upstream rest 
frame is initially a beam-like or ring-like distribution. 
These velocity distributions excit e electromagnetic field 
and am plify the magn etic field (Wu & Davidson' '1972 ; 
iLee fc Id.. 1987; .Ravmond et al.l [2008; .Ohira ct al...2009 : 
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lOhira fc Takaharal [20I0I ) . This is a promising mecha- 
nism to explain s ome observations conce r ning to the large 
magnetic field (iVi nk fc Laming* '2003'; 'Bere zhko et alj 
120031 iBamba et alT 2005; Uchiyama ct al. 200^. More- 
over, the electromagnetic instabilities would heat the up- 
stream region. Ha observed from the upstream region 
(|Lee et alhoWt) can be interpreted as the results of leak- 
age neutral particles. 

In Section 13.11 and 13. 2[ we assumed that the adia- 
batic index of pickup ions is 5/3. However, there is 
no guarantee that the behavior of the pickup ions is 
the same as the standard gas because of the collision- 
less system. Especially, the behavior of pickup ions at 
the shock is important for the shock jump condition and 
particle accelerations. Even though the Mach number at 
the shock is not so small, the compression ratio could be 
smaller than 4 because the adiabatic index of picku p ions 
could be larger than 5/3 f|Fahr fc Chalov 2008; Wu et al.l 
I2OO9I) . If pickup ions drain the large fraction of the 
shock kinetic energy, upstream ions are not heated up 
to T = 3mM L/16. This can ex plain the recent observa- 
tion of Ha jlielder et all 120091 ) which showed that the 
temperature derived from the line width of Ha is much 
smaller than derived from the proper motion. These open 
issues are crucial for not only particle accelerations and 
dissipation but also the amount of leakage neutral parti- 
cles from the downstream region to the upstream region. 
Therefore, we treated values concerning to the leakage 
neutral particl es as free parame ters in this paper. On the 
other hand, Blasi et al.l (j2012l) provided the precise ve- 
locity distribution of neutral particles by assuming that 
only the ion distribution is a Maxwellian. 



In this paper, we did not specify neutral particles. He- 
lium atoms have a smaller cross section than that of hy- 
drogen atoms. We expect large scale of the neutral par- 
ticle precursor compared with that of hydrogen atoms. 
However, the ionization fraction of helium depends on 
time because helium ato ms are ionized by radiat ion from 
the downstream region dGh avamian "eFaL|[2000l) . There- 
fore, the precursor scale and the injection of helium ions 
into DSA could depend on an SNR age. The CR injec- 
tion history of helium ions is important to understand 
the observed CR helium spectrum at the Earth (Drur^ 
IMllOhna fc Iokall20Tll) . 

In this paper, we have investigated the effects of leak- 
age neutral particles on shocks by using four fluid ap- 
proximations (upstream neutral particles, upstream ions, 
leakage neutral particles from the downstream region and 
pickup ions) . We have calculated the precursor structure 
due to the leakage neutral particles and found analyti- 
cal solutions by using the decoupling approximation or 
the tight-coupling approximation, where the decoupling 
means that upstream neutral particles do not interact 
with upstream ions and the tight-coupling means the op- 
posite case. We found that for small leakage of neutral 
particles, the shock compression ratio becomes signifi- 
cantly small. This can explain the observed gamma-ray 
spectra slightly steeper than the simplest prediction of 
DSA. 
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